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1. Introduction 

The top quark was discovered in 1995 and has been intensively studied at the Teva- 
tron. It will also be a prime object of study at the forthcoming Large Hadron Collider 
(LHC). With production rates far in excess of current experiments, the LHC can be 
thought of as a veritable top quark factory. The large mass of the quark, near the 
electroweak scale, enables detailed scrutiny of its interactions unshrouded by hadroni- 
sation effects (which are suppressed by powers of the perturbative momentum scale). 
Furthermore, effects of physics beyond the Standard Model are thought to lie around 
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the electroweak scale, thus the top quark sector provides a valuable window for new 
physics. 

A useful process to study in this regard is the production of single top quarks 
via the weak interaction, as it offers a relatively clean probe of the properties of 
the heavy quark. Although there was recently at least a 3<r evidence of single top 
production at both Tevatron experiments (see e.g. refs. [1, 2]), the cross section is 
sufficiently small (within the SM) as to preclude detailed scrutiny. However, one 
expects a significant number of single top events at the LHC, where the centre of 
mass energy is much higher. 

Accurate estimates of rates and kinematic distributions at hadron colliders ne- 
cessitate higher-order perturbative computations in QCD. Furthermore, in order to 
optimise acceptance cuts in experimental analyses or perform full detector simula- 
tions, one needs realistic hadron-level events. These are obtained using Monte Carlo 
event generators that incorporate the simulation of parton showers and hadronisation 
models. The complementary benefits of fixed-order computations and parton shower 
simulations are by now well-known, as are the advantages of combining them into a 
framework which utilizes the benefits of each of them. The MC@NLO approach [3, 4] 
provides a way of achieving this, by allowing one to match cross sections computed at 
NLO in QCD with an event generator. No modifications to the latter are necessary, 
and colour coherence is preserved by the matching procedure. Therefore, existing 
parton shower Monte Carlos with no add-ons can be used for this purpose. 

There are three distinct production modes for single top quarks. In earlier 
work [5] we included two of these modes, the s- and t-channel processes, into the 
MC@NLO framework. In this paper, we implement the production of a single top 
quark in association with a final state W boson, thereby completing the description 
of single top processes in hadronic collisions at this level of accuracy, including spin 
correlations of decay products, which we incorporate as explained in ref. [6]. 

As is well known, the inclusion of higher order corrections for the Wt channel 
has challenging peculiarities due to interference with the ti process. This interference 
becomes extremely large in certain phase-space regions, and apparently renders the 
perturbative computation of the Wt cross section meaningless. Nevertheless, owing 
to distinct features of ti production, several definitions of the Wt channel have been 
given in the literature, each with the aim of recovering a well-behaved expansion 
in a s . The problem of interference in fact affects any computation that considers 
contributions beyond the leading order, i.e. at least 0{g^ N al). The cross section at 
this order has been previously presented in refs. [7, 8, 9], where only tree-level graphs 
were considered, and in refs. [10, 11, 12], where one-loop contributions were included 
as well. In ref. [13], the calculation of NLO electroweak effects has been carried out, 
using the results of ref. [11] as a basis. We shall comment on some of these papers 
in the following. 

The aim of this paper is to critically examine the definition of the Wt channel, 
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and to propose two options that can be used in the context of an NLO calculation 
interfaced with parton showers. Although these two options are well defined even in 
the absence of kinematic cuts, we will in addition consider final-state cuts that will 
further help the separation of the Wt and ti processes. Such cuts, at variance with 
those of previous approaches (except ref. [9]), are easily applied in an experimental 
environment. The implementation of two definitions of the Wt channel in the same 
framework will allow us to estimate the theoretical systematics potentially affecting 
the extraction of the Wt signal from data. Although the implementation of any 
definition of the Wt channel into any higher-order computation is technically non- 
trivial (and MC@NLO is no exception to that), we argue and make it plausible that 
an NLO+parton shower framework is uniquely suited to the discussion of the physical 
consequences of such a definition. 

The structure of the paper is as follows. In section 2 we give a first general dis- 
cussion of the problem of interference with ti production. In section 3 we recall the 
various modes of single top quark production, before giving some technical details 
on the implementation of Wt production in MC@NLO, and on the underlying NLO 
computation. In section 4 we analyse the problem of interference in more detail, dis- 
cuss the approaches in the literature, and propose two definitions of the Wt channel 
that can be used in a parton shower context and in an experimental analysis. We 
proceed in section 5 to present the results obtained by applying these two defini- 
tions in MC@NLO simulations. We discuss physical implications and conclude in 
section 6. Certain technical details are collected in the appendices. 



2. Nature of the problem 

In the perturbative computation of Wt production, one must consider all possible 
partonic processes with final states 

t + W + ^X t . (2- 1 ) 

i 

Here, {X{\ is a set of particles (partons in a QCD computation), whose multiplicity 
increases as the perturbative order increases. At the leading order (LO) in the SM, 
0(gl / a s ), such a set is empty, and the underlying partonic process is 

bg — > tW . (2.2) 

When next-to-leading order (NLO) corrections in a s are considered, contributions 
(e.g. gg — > tWb) will appear such that 

{Xi} = b . (2.3) 

Some of the relevant Feynman diagrams will feature the Wb pair originating from 
a i internal line; in other words, the momentum flowing in the t propagator will 
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be kw + ki. Therefore, in the computation of certain observables, one will need to 
integrate over the region 

M? w = (k w + k- b ) 2 ^ml (2.4) 

When this is the divergence is encountered, which is regulated only by using 

a finite width for the (anti-)top, T t ^ 0. However, a non-zero top width is an all- 
order result in g w . Thus, the inclusion of higher-order QCD corrections forces one 
to include electroweak corrections to all orders so as to avoid divergences, and this 
potentially spoils the power counting in g w , according to which eq. (2.2) is the LO 
contribution to tW production. 

It turns out that indeed the region of eq. (2.4) causes severe problems, and as 
a result NLO QCD corrections are much larger than the LO result obtained with 
eq. (2.2). One possible way out is that of considering only observables which are 
exclusive enough to allow one to exclude, through final-state cuts, the resonance 
region eq. (2.4). This is easily done in a fixed-order, parton-level theoretical com- 
putation, where the W and the b are easily accessible. It becomes however a more 
indirect procedure in the context of a parton shower simulation, and impossible in a 
real experiment. Note also that, in order to avoid biases in an MC simulation, very 
loose cuts or no cuts at all will have to be imposed in the computation of underlying 
matrix elements, with the result that the majority of the events generated, being 
close to the t resonance, will actually be thrown away by final-state cuts, leading to 
a very low efficiency. 

The dominance of the t resonance suggests that a possible approach is that 
of simply considering W + W~bb final states (with the possibility of also including 
production spin correlations as well, for the fermions resulting from the decays of 
the ly's - see e.g. ref. [14] for di-leptonic decays). In this way, what has previously 
been denoted as the LO contribution to Wt production, eq. (2.2), can be seen as 
a correction to W + W~b observables (i.e., observables inclusive in the final-state 
b present in the W + W~bb matrix elements). The way in which this contribution 
is taken into account is a matter of careful definition, which is especially delicate 
(because of the double counting problem) when interfacing the matrix elements to 
parton showers - see what was done in ref. [9]. 

By emphasising the role of W + W~b{b) final states, these approaches have an 
immediate connection with data, and finding the Wt "signal" is thus a matter of 
a careful counting experiment. The key question is: are these predictions accurate 
enough for the counting to be reliable? Indeed, in these approaches the problem of 
full NLO corrections, whose knowledge is known to be crucial in top physics, is not 
considered. 

The computation of NLO corrections to Wt production in QCD can only be se- 
riously undertaken by recovering a meaningful definition of a perturbative expansion 
whose LO contribution is that of eq. (2.2). We have seen that the problem arises 
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(1) (2) (3) 

Figure 1: Leading order diagrams for single-t production in the (1) s-channel, (2) i-channel 
and (3) Wi-channel. The t-quark line is doubled. 

because of the interplay between a s and g w expansions, and specifically because of 
the necessity of considering all-order contributions in g w . The key observation is that 
the all-order result in g w that we need is associated with a decay, whereas we are only 
interested in the role of EW interactions in the production process. A meaningful 
expansion in a s of the Wt cross section could therefore be achieved if electroweak 
effects in production and decay could be disentangled. Such a procedure can never 
be fully consistent theoretically, but it can be given an operational meaning, which 
can be tested experimentally. The rest of the paper will be devoted to separating 
production and decay EW effects in the context of MC@NLO. 

3. Wt production at NLO 

In this section we introduce the various single top production modes at LO, before 
giving some details of the calculation of the Wt channel at NLO in QCD. Technical 
details regarding the implementation in MC@NLO are presented, although as these 
are similar to those relevant to the other single top production channels [5] we refer 
the reader to previous publications where appropriate. Since in this section we will 
be mainly concerned with the computation of the matrix elements entering the NLO 
predictions, and with the construction of the MC subtraction terms for MC@NLO, 
we can safely ignore the issues discussed in sect. 2, but will return to them in sect. 4. 

3.1 Born level 

There are three modes for the production of a single top quark at LO 1 . Each of 
these provides a separate and complementary means of studying the Wtb vertex: 
the s-channel mode (fig. 1(1)) involves a time-like, off-shell vector boson which may 
reveal new resonances; the t-channel mode (fig. 1(2)) involves a spacelike boson and 
is sensitive to flavour changing neutral currents (FCNC's). Finally, the Wt mode 
(fig. 1(3)) has an on-shell W boson, and thus provides a complementary source of 

1 We assume throughout the paper that a 5 flavour scheme can be used for the quark sector, 
where the b quark is included in the initial state parton distributions. In principle the calculation 
could also be formulated in a 4 flavour scheme where all b quarks are generated in the final state. 
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information on the flavour structure of the vertex (i.e. Vtb), and its chiral properties. 
Based on expectations from the SM, this latter process has far too small a cross 
section to be observed at the Tevatron, but is expected to be significant at the LHC. 
As shown in fig. 1(3), the lowest order partonic process for Wt production is 

6(pi) + 9{P2) — > t(fci) + W~(k 2 ) . (3.1) 

Our description is confined to W~t production. As explained in appendix A, for 
W + i production no new diagrams need be computed, so that the description below 
suffices. Where no confusion is possible in the following, we denote the light quark 
attached to the Wt vertex as a b, implicitly representing any CKM-allowed down-type 
quark. 

It is convenient to introduce the following invariants 

■s = (pi+P2) 2 , h = t — mf = (ki-pi) 2 -m 2 , u\ = u-m 2 t = (k 2 -pi) 2 -m 2 t (3.2) 

such that s + ti + Ui — — mf , with m t denoting the top quark mass. We treat 
all other quarks as massless, including the b quark. The lowest order cross section 
can be written 

d<r<® = M^dh , (3.3) 
M {0) = — - - q 2 ^ U (0) l 2 (3 4) 

where g s , g w are the QCD and EW coupling constants; the number of colours; 
dfo denotes the two-body final state phase space, and the spin-summed Born-level 
matrix element is given by 

+ 2 m 2 w (s + u^Uj + m 2 t m 2 w s ( _ m 2 _ m\ \ 
su\ V 2m w 2m w) 

sui V m w mw) 

3.2 NLO computation 

In order to implement a process in MC@NLO, one must cancel the singularities 
arising in the real and virtual graphs in the particular subtraction formalism of 
refs. [15, 16] and denoted FKS henceforth. In the FKS formalism, the NLO cross 
section is expressed in terms of the finite quantities that result from the cancellation of 
the (universal and process-independent) soft and collinear poles arising from virtual 
and real corrections. In order to obtain these finite quantities for a given process, one 



-4)- (3-5) 
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needs the real matrix elements computed in 4 dimensions, and the virtual corrections 
computed in d dimensions. Although not strictly necessary in the FKS method, a 
rf-dimensional virtual computation was carried out to fix the convention for its finite 
part. 

The real and virtual matrix elements relevant to Wt production can be extracted 
from the NLO calculation of the Wc process [17], for which the diagrams are almost 
identical 2 . Given that the computation of ref. [17] adopted the phase-space slicing 
method [18, 19, 20] for cancelling the real and virtual singularities, we converted the 
conventions of the phase-space slicing into those of FKS. In order to check that this 
was carried out correctly, we repeated the calculation of the matrix elements. 

3.2.1 Virtual corrections 

As stated above, we calculated the one-loop virtual diagrams in dimensional regu- 
larisation in d = 4 — 2e dimensions. Tensor and vector integrals were reduced to 
scalar integrals using the standard Passarino-Veltman algorithm [21]. For generation 
and evaluation of the relevant amplitudes, we utilised the Mathematica packages 
FeynArts [22] and FeynCalc [23], together with the FORM [24] computer program. 
As discussed previously, the results were checked against the virtual contributions 
obtained for Wc production in ref. [17], and found to be in agreement after the nec- 
essary analytical continuation (from m c < my/ to m t > mw)- We have also checked, 
analytically and numerically, that our results for the scalar integrals are in agreement 
with those of ref. [25]. 

The virtual corrections are a Laurent series in the parameter e, with double and 
single poles arising from soft, initial-state collinear, and ultraviolet (UV) singular- 
ities 3 . To remove UV poles we renormalise the top quark mass using an on-shell 
condition, with the QCD coupling in the MS scheme modified such that the top 
quark loop contribution is subtracted on-shell. This particular scheme [26] ensures 
that the top quark virtual contributions decouple in the limit of small external mo- 
menta. Specifically, the coupling is renormalised through next-to-leading order as 

^ 9s «) [l + (± + - >n 4.) ( A) \ + S£> | , Q)] 

(3.6) 

where /x F is the factorisation scale, and /i R the renormalisation scale. Furthermore, 
Po = {HC a — 2n/)/3 with rif equal to the number of light flavors (here five) plus 
one. From this condition one can derive the following relation for the renormalised 
QCD coupling 

^ = _ M f _M> ( A + 1) + 0(ai) , ( , 7) 

2 Essentially, only the W + (cc) channel has no analogue in the present case. 
3 Final-state collinear singularities are absent due to regularisation by the top mass. 
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which indeed removes the top quark loop from the /^-function, such that the number 
of light flavours implemented in the code is rif. The renormalisation of the top quark 
mass is given by 



m t — > m t + 8m t 
m t 



1 + -^ Cf 



-3 



+ 37 B -31n47r-4-31n 



/4 

2 

mi 



(3.1 



In contrast to the s and t single top production channels, the top quark occurs as an 
internal line in the Born amplitude. This gives rise to a contribution to the amplitude 
coming from the expansion of the renormalised top quark propagator 

i i i „ 1 



y-K + fcn.) t- mt + J^ t Sm 'J^t (3 ' 9) 

The first term in this expansion gives the Born amplitude ^4°, whereas the second 
gives a modified amplitude A'^ involving the "squared" top quark propagator. In 
terms of the finite remainder of 5m t , the cross section receives the contribution 



da 



(hV) 
Smt 



(Sm t ) 



1 1 



1 



finite 



2sAN(N 2 - 1 



-o 2 — 

r s 8 



^(o U ,t(D + A , 



(3.10) 



where 

A (o) A m + A ^) A m = _ 16NCf 



9 * A \4m 6 t s + Amts(rriL + m)] 



+ m\ [-8m^s + 2m 2 v sui + U\{u\ + t 1 )(s + 2ti)] 
+ «i [«iti(ui +h) + m 2 v {-s{u 1 -2ti) +4*i(ui + ti))] }>. (3.11) 



Following the renormalisation procedure, soft and collinear singularites remain and 
are proportional to the Born cross section such that the virtual terms may be written 
(c.f. eqs. (3.1) and (3.2) in ref. [15]) 



daW = C e g 



~(C F + C A ) + -C A [hi^ 
e 2 e \ mf 



In 



-Ml 



• -( V (a In f ' 



5 



mi 



e 



In 



mt 



daf2 2e + da^J, (3.12) 



where da^} 2 e * s the Born cross section in 4 — 2e dimensions, and 



a = 



16tt 2 



m; 



(3.13) 



These remaining poles are cancelled by similar singularities in the real contributions 
and by the collinear counterterms that arise from the renormalisation of the parton 
densities, in a form prescribed by the FKS formalism. 
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3.2.2 Real corrections 

Since only initial-state collinear singularities are present in the process considered 
here, it is not necessary to partition the phase-space according to the FKS prescrip- 
tion, as done in ref. [5] for the case of s- and t-channel single-top production (see 
sect. 2.1.2 of that paper). 

Therefore, the only technical difficulty in assembling the real corrections is in the 
explicit calculation of the diagrams 4 . As explain previously, we used the results of 
the Wc calculation presented in ref. [17]. There, the matrix elements were computed 
by considering the fictitious W boson decays 

W(q) — Q(p) + b 3 (q 3 ) + g A (q A ) + g 5 (q 5 ) (3.14) 

and 

W(q) — Q(p) + b 3 (q 3 ) + b 4 (q 4 ) + b 5 (q 5 ) , (3.15) 

where Q denotes generically the only quark with mass different from zero. These were 
decomposed into colour-ordered amplitudes, and computed using FORM [24]. The 
Wc and Wt matrix elements were subsequently obtained by crossing and summing 
over colour orders. This saved somewhat on computational effort, since the colour- 
ordered amplitudes served a second purpose, namely in the accounting of colour 
connections in the parton shower stage of the MC@NLO construction. 

As a further check on the real computation, we compared the results we obtained 
from our calculation against the corresponding tree level matrix elements generated 
using the MADGRAPH program [27, 28]. We found agreement in all cases. 

3.3 Implementation in MC@NLO 

Several processes have by now been computed in the MC@NLO formalism, and we 
will therefore refrain from describing here the necessary steps for the implementation 
of a new reaction. The interested reader can find all details in refs. [3, 4, 5], which 
also report all the relevant analytic formulae we need (these are process-independent, 
and therefore no new computation is specifically required for Wt production). In this 
section, we will limit ourselves to giving the necessary information for the construc- 
tion of the MC subtraction terms (which are process dependent) necessary for the 
matching of the NLO computation with HERWIG [29, 30]. 

One may write the MC subtraction terms as (see refs. [4, 5]) 

(3.16) 



4 There is a subtlety involving the construction of the local initial-state collinear counterterms 
used in numerical codes, as in all cases in which a gluon is exchanged. This issue is discussed in 
appendix B. 
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, -EEE^f" 
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where % sums over the different partonic subprocesses, and L e {+, — , /i} labels the 
parton leg from which the FKS parton is emitted (respectively: incoming parton 
along the +z direction; incoming parton along the — z direction; final state top 
quark). The index / runs over the different colour structures. If q a and q@ denote 
the 4-momenta of the colour partners relevant to a given branching, the shower scale 
associated with the branching is 



El 



\q a -qp\- 



(3.17) 



The individual subtraction terms on the right-hand side of eq. (3.16) have the form 
(see eqs. (5.2)-(5.5) in ref. [4]) 



da. 



(+,0 



da. 



(-.0 



MC % 



dx u dx 2i , 



MC 



dxu dxu i 



da 



(h,i) 



f^if)A H2 \^f)da^ 



dx±f dx2f , 



(3.18) 
(3.19) 
(3.20) 



where f^ 1 ' 2 ^ are the initial state parton distributions, and zl labels the partonic 
momentum fraction carried by the FKS parton. The short- distance cross sections 
da[ L ' 1 ^ can be obtained from eqs. (5.6)-(5.8) of ref. [4] 



da. 



(±>0 



da. 



= gS^,(4?)^e(i-ef)ef4?- 

Vi V 




where M = 1 for quark branchings, and 2 for gluon branchings (since the gluon has 
two colour partners); see refs. [4, 5] for the definition of da. The various partonic 
subprocesses contributing to the Wt channel are listed in table 1 (the initial states 
bq and qb represent also bq and qb). The shower scales to be used in eqs. (3.21) 
and (3.22) are equal to the absolute values of the dot products given in table 1. 
From table 1, we see that Af = 2 for all branchings in the processes with {b, q}, {b : b} 
and {b, b} initial states, and for branchings from leg — and leg + in the processes 
(6, g; t, W, g) and (g, b; t, W, g) respectively. 

When interfacing the NLO calculation with a parton shower, it is also necessary 
to supply a colour flow to the Monte Carlo. When a given subprocess can have more 
than one colour flow, the one that is given to the parton shower is chosen on a statisti- 
cal basis. The weights that govern these probabilities are formed using leading order 
(in the number of colours N) approximations to the squared matrix elements of the 
corresponding tree-level diagrams. In our case, these have already been calculated, 
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NLO 


{o, g, 1, vv ) 


in h- t W\ 

\9i °, r , vv ) 


(9,9 


1 T I 7" 7 \ 

t, W, 6) 


+iprp2) 


{ — - \ 
-(P1-P2) 


(6,5; 


t, W, g) 


/ — — — T \ 

-{PrP2,P2-h) 




(5,6; 


J. T T T \ 

t, W, q) 




— 

+ (prP2,Prh) 


(6,5; 


t,W~b) 


_ , 

-(Pl-P2,P2-h) 




(5,6; 


t,W,b) 




+ (prP2,Prh) 


(6,6; 


t, W, b) 


-(pl-P2,P2-h) 


+(pi-p2,Pi-h) 


(b,g 


t,W,g) 


+, fi(prh); -(pi-p2,P2-k 1 ) 




(9,b 


t,W,g) 




-,fi(p2-h); +(prP2,Prh) 



Table 1: Short-distance contributions to MC subtraction terms. The two columns corre- 
spond to the two possible Born cross sections. For a given process, the entries show the 
emitting legs, and in round brackets the value(s) of the shower scale(s) Eq (up to a sign). 



as the real emission contributions were computed using colour ordered amplitudes. 
Given that there is a one-to-one mapping between the colour ordered amplitudes and 
the various colour flows, one can choose to use these colour ordered amplitudes in 
the statistical determination of the colour flow 5 . 

4. Interference between ti and Wt production 

In section 3 we have discussed the calculation of the Wt cross section within the FKS 
subtraction formalism, and its subsequent implementation in MC@NLO. However, 
the use of such a Monte Carlo event generator assumes that the Wt channel is well- 
defined beyond LO. In fact, as already outlined in sect. 2, this is not the case, and 
the theoretical definition of this production channel is not straightforward. 

At LO, the Wt pair is produced through the reaction of eq. (3.1), and its cross 
section is smaller than that of top pair production by a factor of about 15. Beyond 
LO, some of the Feynman graphs that contribute to the Wt channel and that are 
dominant in the region of eq. (2.4) are shown in fig. 2. These diagrams can be 
interpreted as the production of a tt pair at LO, with subsequent decay of the t into 
a bW pair. It follows that the set of Feynman graphs contributing to gg — > tWb or 
qq — > tWb can be divided into two subsets, that are customarily denoted as doubly 
resonant and singly resonant, the former being those depicted in fig. 2 or analogous 
ones. The problem discussed in sect. 2 can be formulated in terms of the interference 
between these two subsets of graphs, and physically interpreted as the interference 
between Wt and tt production. This interpretation is consistent with our aim of 

5 Other choices are possible, as one is free to modify the above procedure by terms sublcading in 
the number of colours N. 
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Figure 2: Diagrams that are doubly-resonant, in the sense that the intermediate t can be 
on-shell. 

emphasising the role of electroweak interactions in the production mechanism of the 
t and of the W . The interference would not cause any problem if the contribution 
from the doubly-resonant diagrams was small. However, a large increase of the cross 
section occurs when the t propagator becomes resonant i.e. in the kinematic region 
of eq. (2.4). It follows that the naive calculation of the Wt channel, defined as that 
resulting from the straightforward computation of Feynman diagram as discussed in 
sect. 3, becomes meaningless in perturbation theory. In this section, we will discuss 
how it is possible to separate out Wt and ti production, and recover a workable 
definition of the Wt channel. 

4.1 Previous approaches 

This interference problem has already been discussed in several papers. Given that 
the problem here only arises for gg and qq initial states, present in real-emission 
processes, it is not peculiar to full higher-order computations, but also occurs in 
those that only include tree-level graphs. Each of the previous Wt computations 
beyond LO has introduced a strategy for isolating the Wt production channel. We 
briefly discuss each of them in turn. 

In ref. [8], the suggestion was made to place a cut on the invariant mass M^ w of 
the Wb pair, which can be written in the form: 



where r t is the width of the top quark. This cut thus removes events from the re- 
gion of phase space corresponding to the t resonance. It efficiently reduces the ti 
interference to the Wt channel, but cannot be directly applied in an event gener- 
ation/experimental situation, where one is unable to unambiguously identify the b 
quark and W boson which originate in the hard interaction. 

In ref. [7], the Wt signal is defined by subtracting, at the level of squared am- 
plitudes, the ti cross section multiplied by t — > Wb branching ratio; the procedure 
is defined in a fully-inclusive way. A comparison was made with the invariant mass 
cut approach of ref. [8], and the resulting cross section found to agree if the choice 
k ~ 15 was made. 



M- bW -m t \> kT; 



(4.1) 



12 



As discussed in sect. 2, the approach of ref. [9] addresses the problem in a more 
general context, i.e. that of simulating W + W~b(b) production in a realistic parton- 
shower environment. As such, the definition of the Wt channel does not need to 
be directly considered, except for the necessity of avoiding double counting. This 
way of by-passing the interference problem is correct, and facilitates the comparison 
with data. However, it does not allow one to include the full NLO corrections into 
the results, because of both problems of practice (the relevant one-loop, six-point 
amplitudes - or eight-point amplitudes, if fully decayed W's are considered - are 
not available yet), and of principle (how to combine the one-loop contribution for 
W + W~b with the rest). 

The full NLO computation of ref. [10] follows the same type of strategy as ref. [7], 
where the total cross section is modified by a subtraction term which effectively 
removes the tt-like contribution. However, few details on how this subtraction term 
is defined are given. The results are presented only for total rates. 

Finally, a definition of the Wt channel within a fully differential NLO compu- 
tation was presented in ref. [11]. There, the suggestion was made to place a cut on 
the transverse momentum p T of the "additional b quark that appears at the next-to- 
leading order" [11], namely the b quark that appears in the diagrams of fig. 2, which 
accompanies the b quark coming from the decay of the top (not shown in fig. 2). In 
practice, this amounts to keeping events that satisfy 

p\<p { r°\ (4.2) 

where appropriate. This means that in the case of partonic processes with no such 
additional b quark, no veto is applied. The reasoning for the condition in eq. (4.2) 
is that harder b quarks tend to have come from the decay of a top, and thus the 
probability for producing events with two hard b quarks is dominated by the contri- 
bution of diagrams with a t and a t both on-shell or almost on-shell. By requiring 
the additional b to be softer, one thus reduces the interference between the doubly 
resonant diagrams and the singly resonant diagrams which are identified with the 
Wt channel. Furthermore in ref. [11], the factorisation scale is chosen to be equal to 
p(veto) ^ j£ condition is not met, doubly-resonant diagrams are removed from the 
computation at the amplitude level. Finally, processes with a qq initial state are not 
included, independently of the choice of factorisation scale. This is only a problem 
of principle, since in practice the numerical impact of such processes is so small that 
they can be neglected. As in the other calculations discussed above, it is not possi- 
ble to apply this definition of the Wt channel in an event generation/experimental 
situation, since the veto is based on a partonic picture at a given order in a s . 

4.2 MC@NLO approach 

In this section, we consider the problem of defining the Wt channel in a way that is 
applicable in an event generator context, where both initial- and final-state parton 
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showers are present. Although none of the approaches discussed in subsection 4.1 can 
be directly applied in this situation, they can be used to some degree in motivating 
suitable solutions. 

The simplest and most drastic solution is to remove from the computation the 
contributions of those processes which contain doubly-resonant diagrams. At the 
NLO, this amounts to removing completely diagrams having gg or qq initial states, i.e. 
regardless of whether or not they are doubly resonant. As far as the qq contribution is 
concerned, this is what is done in MCFM [11]. It is not quite the same as what is done 
for the gg initial state, but is numerically very close to that, owing to the combination 
of the veto, and of the choice of the factorisation scale /x F = p^ eto ^ . However, even if 
one generalises the definition of the veto to apply in a parton shower context (as we 
will do in the following), the removal of processes characterised by a particular initial 
state is theoretically unfeasible. Firstly, renormalisation group invariance is violated 
such that the underlying NLO calculation has the same formal scale dependence as a 
LO one (thus undermining one of the motivations for using an NLO generator at all). 
Secondly, and perhaps more seriously, partonic processes mix as soon as higher-order 
effects are considered, as happens when initial-state showers are present. Thus, we 
discount this method as a viable means of separating out the Wt channel in a full 
Monte Carlo generator and do not present further results from it here. 

Instead, we present two definitions of the Wt channel, that are designed in such 
a way that, by comparing them, one can directly assess the impact of the interference 
with ti. Thus, if the results from the two definitions agree, we can be confident of 
having isolated the Wt channel. Such an agreement, when it occurs, is the result 
both of the definitions themselves, and of final-state cuts which may or may not be 
applied. It is important to stress that the definitions that we shall give are meaningful 
even without any subsequent cuts. This engenders a greater degree of flexibility in 
the practical investigation of whether the Wt channel is well defined or not. 

Our two definitions can be summarised as follows: 

1. Diagram Removal (DR). Here one simply removes all diagrams in the NLO Wt 
amplitudes that are doubly resonant (i.e. those diagrams shown in figure 2). 

2. Diagram Subtraction (DS). In this approach, one modifies the NLO Wt cross 
section by implementing a subtraction term designed to cancel locally the ti 
contribution. 

Note that DR is different from the removal of gg- and gg-initiated processes discussed 
previously, since diagrams with these initial states are kept if they are not doubly 
resonant. Note also that the DS procedure is similar to what has been proposed 
in ref. [7], but the construction of its local subtraction term involves some technical 
complications, which we discuss later. 
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In order to discuss in detail the DR and DS definitions, we introduce some 
notation. Let us denote by A a /3 the 0(g w a s ) amplitude of the process: 

«(pi) + Pfa) — > t(h) + W(k 2 ) + 5(k) . (4.3) 

The precise identity of parton S is not relevant for what follows, so we omit it from 
the notation. We have: 

where the two terms on the r.h.s. are the respective contributions of the singly- and 
doubly-resonant diagrams to the naive Wt cross section. In the computation of the 
cross section at NLO, the following quantity appears: 



1 2 AWt) 



+ 28{ylS'U<'f} + |4?| (4.5) 

<Saf3 + T a p + V a p . (4.6) 



We stress that the terms T a p and T> a p are non-zero only in the cases of {a, f3} = {g, g} 
and {q,q}, where in the latter process q may also be a b quark. In these cases, 5 
will be a b quark or another down-type antiquark according to the CKM matrix. 
Furthermore, V a p has neither soft nor collinear singularities, while those of T a p are 
integrable and subleading w.r.t. those of S a p. Therefore, in the context of an NLO 
computation in the FKS formalism (or for that matter in any subtraction formalism), 
the NLO real-emission contribution to the subtracted short- distance partonic cross 
section including the flux factor and phase space will be: 

dcr a p = y s {paP + J «/3 + V ^p) dfo , (4-7) 

where the hat denotes that infrared singularities have been suitably subtracted. It 
is understood that the subtraction is performed by means of plus-type distributions, 
which therefore may apply to the phase space as well. The hadroproduction cross 
section resulting from eq. (4.7) is 

da = dcr^ + / dx 1 dx 2 C a /3da a/3 

a/3 J 

= daW + |^§£"/3 (4/9 + T aP + V a ^j # 3 , (4.8) 

a/3 

with S the squared centre of mass energy of the colliding hadrons and £ Q/ g the parton- 
level luminosity. The quantity da^ 2 \ requiring a two-body phase space, denotes all 
contributions to the cross sections that are not already included in eq. (4.7), i.e. the 
Born, soft-virtual, and collinear remainder terms. Both da^ and eq. (4.7) (but not 
their sum, eq. (4.8)) are convention-dependent, since finite pieces can be freely moved 
from one contribution to the other, but this is irrelevant for the following discussion. 
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When the NLO computation is then matched to parton showers according to the 
MC@NLO prescription, the above equation must be modified by the subtraction of 
MC counterterms. We can choose to absorb these in S a p, because this is the only 
piece that contains leading soft and collinear singularities. Thus the schematic form 
of eq. (4.8) applies at both the NLO and MC@NLO levels. In this notation, the DR 
cross section corresponds to: 

<*7< D *> = da^ + W t^C^dfo , (4.9) 

a/3 

i.e. there are now no terms T a p or T> a/ 3, as all doubly resonant diagrams have been 
removed from the amplitude. As mentioned previously, this cross section violates 
gauge invariance; this issue will be discussed in sect. 5.2. 

Starting from eq. (4.8), we also define the DS cross section. This amounts to 
writing: 

da {us) = da - da sub \ (4.10) 

where da subt is designed to remove numerically the doubly-resonant contribution. 
This may be achieved locally by defining 

da suu = J2j dx,dx 2 C a/3 da s a f ; (4.11) 

a/3 

da s a f = j-V a pd<p 3 , (4.12) 

such that the quantity 

V a(3 -V af3 (4.13) 

will vanish when M^ w = (k + k 2 ) 2 mf . Note that V a/3 and V a p themselves will, 
in such a limit, either diverge, if r t = 0, or have a Breit-Wigner-like peak, if r t ^ 0. 
The DS cross section in eq. (4.10) can now be re- written in the same form as eq. (4.9): 

da^ = da® + W |r|| £ ^ + ^ + V aP - Vj) dfo . (4.14) 

a/3 

One sees that the difference between the DR and DS cross sections has the form: 

da^ - da<r*> = j + V*fi ~ dfo , (4.15) 

a{3 

and thus is composed of a contribution from the interference term, and of the differ- 
ence between the subtraction term and the true doubly resonant contribution to the 
NLO cross section. 

Our aim is now to construct a gauge-invariant subtraction term, such that the 
difference V a/ 3 — T> a p is as close to zero as possible. Note also that requiring the 
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subtraction term to be local and gauge invariant prevents this difference from being 
identically zero. The subtraction term should have the schematic form: 



V a/3 = \A (0 \a(3 ^tt)\ 2 x BW(M- bW ) x \A {0) (t -> Wb)\ 2 (4.16) 

where A^\a(3 — > ££) is the LO (i.e., 0(a s )) top pair production amplitude, A^°\t — > 
W&) the decay amplitude of the anti-top, and BW(Mg W /) is the Breit-Wigner func- 
tion. 

The first obvious difficulty in constructing a workable implementation of eq. (4.16) 
is that the kinematics on the l.h.s. should be that of the full a (3 — > tWb process (due 
to the requirement of local cancellation in eq. (4.14)), whereas at the same time the 
t needs to be on-shell in order to compute the quantity A^\af3 — > tt) in a gauge- 
invariant way. A second difficulty is that, while A^\a(3 — > tt) is computed with top 
width T t = 0, the Breit-Wigner factor in eq. (4.16) requires non-zero r t . 

As a first attempt to overcome these difficulties one can reshuffle the momenta 
of the decay products, to obtain an on-shell t quark. In this case, one uses a non-zero 
top width only in the Breit-Wigner factor. However, a more fundamental problem 
with eq. (4.16) is that spin correlations are not included, thus spoiling again the local 
cancellation property we seek to achieve. Therefore, rather than eq. (4.16), we need 
a subtraction term of the form: 



— 



A m 



(4.17) 

reshuffled 



where A^j is defined in eq. (4.4). In other words, V a/3 would be identical to T> a/3 , 
were it not for the reshuffling. Being a full amplitude, eq. (4.17) does implement spin 
correlations in the decay of the t. A divergence present in eq. (4.17) when r t = 
disappears by setting r t ^ 0. However, the reshuffling implies that V af3 is not likely 
to have the Breit-Wigner shape in M bW that would be desirable in order for the 
difference in eq. (4.13) to be as small as possible in the whole phase space. This is 
easily rectified by defining: 



BW(M- b 



a(3 ~ BW{m t ) 



2 



(4.18) 

reshuffled 



By construction, the amplitude of this subtraction term at M bW = m t is precisely 
such as to cancel the resonant contribution to the NLO Wt cross-section. 

However, even after this modification, there is a problem. For appropriate cancel- 
lations of collinear singularities, we would also have to use T t ^ in the computation 

l 2 

of \Aai3\ ■ Unfortunately, after setting T t ^ in the radiative amplitudes in this way, 
their collinear limits are modified, and in order not to disrupt the local cancellation 
of collinear divergences one would need to perform the full calculation adopting a 
framework for the consistent inclusion of finite-width effects in all kinematic regions 
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of interest (e.g. the complex mass scheme [31]). Given the scope of the present pa- 
per, a more pragmatic approach suffices, which has been extensively used at LEP in 
dealing with resonant decays. Namely, we set T t ^ only in the doubly-resonant dia- 
grams, and leave Y t = in the singly-resonant ones. This strategy is straightforward 
to implement in MADGRAPH. 



5. Results 

In the previous section we discussed the problem of interference between tt-like am- 
plitudes (followed by decay of the t) and I¥t-like amplitudes. Two separation mech- 
anisms, DR and DS, were defined as being suitable for an all orders computation 
of the scattering amplitude, as occurs in a Monte Carlo event generator. We have 
implemented both of these separation mechanisms in the MC@NLO framework, as 
discussed in sect. 3. We stress that the MC subtraction terms given in sect. 3 are 
identical for DS and DR, since such terms only modify (w.r.t. a pure-NLO compu- 
tation) the form of <S Q( g, which is identical in the DS and DR cross sections. In this 
section we present sample results, and compare in detail the output of the DS and 
DR calculations. Our aim is to analyse, from a perturbative point of view, the degree 
to which separation of the Wt channel is possible, thus giving an upper bound for 
the impact of the interference that can be obtained with a realistic analysis. De- 
tailed phenomenological results will not be presented here, but are postponed to a 
forthcoming publication. In particular, we only consider fully leptonic decays of the 
final-state W bosons, as these are sufficient to furnish a comparison between the 
different separation mechanisms. In practice, semileptonic rather than fully leptonic 
decays will be studied first by experiments at the LHC. Furthermore, our results are 
obtained by neglecting production spin correlations, which is not restrictive as far as 
comparing the two definitions of the Wt channel given here is concerned. We have 
however implemented spin correlations in the DR calculation, using the method of 
ref. [6]. Their implementation in the DS calculation along those lines is technically 
slightly more complicated (but possible), owing to the negativity of the "squared" 
matrix elements at some phase-space points (due to the subtraction term), and is 
also deferred to the future. 

All of the following numerical results have been obtained for the LHC with the 
MRST2002 default PDF set [32], setting the top mass and width to m t = 170.9 
GeV and T t = 1.41 GeV, as well as the W mass and width to m\y = 80.4 GeV and 
T\y = 2.141 GeV. The default values of the renormalisation and factorisation scales 
are equal to the top mass. The LO results quoted in this section have been obtained 
using the same parameters (including PDFs and two- loop a s ) as those adopted for 
NLO computations. 
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5.1 Transverse momentum veto 



As discussed previously, the definitions we gave of the Wt channel are independent 
of any cuts that are subsequently used to further reduce the interference from tt 
production. Nevertheless, it is unrealistic to define the Wt channel, in an actual 
experimental environment, with no cuts 6 at all. It is not the purpose of this paper to 
undertake a thorough investigation of various cut strategies. Therefore, we shall limit 
our considerations to one particular cut, motivated by the transverse momentum veto 
proposed in ref. [11], and discussed here in sect. 4.1. As stated there, the veto acts 
on the additional b quark that may be present at the NLO level. Here we generalise 
this idea to an event generator context i.e. to a situation in which we cannot tell 
with certainty which parton or hadron is associated with the additional b quark. 

Firstly, one searches all ^-flavoured hadrons in a given event, and orders them in 
transverse momentum (or transverse energy). All hadrons whose pseudorapidity t]b 
is outside a given range, which we choose to be 



are ignored. Secondly, a veto analogous to that of eq. (4.2) is applied on the transverse 
momentum (or transverse energy) of the second hardest B hadron satisfying eq. (5.1): 



In cases where no second-hardest B hadron satisfying the above requirements can be 
found, the event is accepted. Note that this is the case in processes which do interfere 
with tt production, but in which a b quark is replaced by another down-type quark - 
an effect off-diagonal in the CKM matrix. In this paper, we assume a 100% 6-tagging 
efficiency However, it is clear that the veto procedure proposed here can be applied 
in more realistic 6-tagging scenarios. A more detailed phenomenological investigation 
of this issue is postponed to a future publication. 

In a realistic analysis the veto would be accompanied by a number of further cuts. 
For example, in the semileptonic decay mode, the veto can be used in combination 
with a jet topology cut, to further reduce the contamination of the Wt "signal" from 
tt "background" . Given that such cuts are absent in the results presented here, our 
findings correspond to a somewhat pessimistic scenario as far as the purity of the 
Wt signal is concerned. 

The definition of the veto adopted here can also be used in the context of a 
pure-NLO, parton-level computation, such as that of ref. [11], by simply replacing B 
hadrons with b quarks 7 . It is instructive to see, however, that the veto is a much more 

6 We mean here cuts whose only purpose is that of separating the Wt and ti processes, and not 
generic selection cuts aimed at selecting the t and/or the W in a collision event. 

7 b quarks can also be considered in an MC context. We have found that MCQNLO results 
obtained by imposing b quark vetos are very similar to those obtained with B hadron vetos. 



Vb\ < 2.5, 



(5.1) 




(veto) 



(5.2) 
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natural constraint in an event generator context than in a parton-level, fixed-order 
computation. In table 2 we present the results for the DR and DS total cross sections, 



MC@NLO 


NLO 


<j(DR) 






















p( veto) = 1Q 


34.66 


33.89 


26.60 


1.30 


1.27 


35.05 


34.74 


34.67 


1.01 


1.00 




p^ = 30 GeV 


41.86 


40.74 


31.85 


1.31 


1.28 


39.93 


39.67 


34.67 


1.15 


1.14 




p^ = 50 GeV 


44.61 


42.92 


33.71 


1.32 


1.27 


42.81 


42.00 


34.67 


1.23 


1.21 


p[ veto) = 7Q GeV 


45.63 


43.65 


34.31 


1.33 


1.27 


44.41 


42.90 


34.67 


1.28 


1.24 


Pt 


= oo 


46.33 


44.12 


34.67 


1.34 


1.27 


46.33 


44.12 


34.67 


1.34 


1.27 



Table 2: Results for the total DR, DS, and leading order cross sections (in pb), obtained 
with MC@NLO (five left columns) and our pure-NLO parton level computation (five right 
columns). See the text for details. The notation p^ eto ^ = oo denotes no veto at all. 



obtained with MC@NLO and with our pure-NLO parton level computation. We also 
give the results (denoted by a^ ^) obtained in the two frameworks by keeping only 
LO matrix elements. In the case of MC@NLO, a^ ^ is thus equal to what one would 
get by simply running HERWIG standalone 8 . The table finally reports the values of 
the ratios of NLO over LO cross sections: 

O-(LO) ' (j(LO) V ! 

In the fixed-order computation at LO, there is simply no second-hardest b quark, and 
o-( L °) is therefore independent of the value of p^ eto ^ . On the other hand, in an MC 
context the initial-state b quark entering the hard partonic process at LO eventually 
results, because of parton showers, in the generation of a B hadron. Therefore, a^ ^ 
computed with an event generator does depend on p^ et0 ^ . The interesting thing about 
this dependence is that it appears to be the same as that obtained at the NLO, as 
can be inferred from the basically constant values of i^( DR ) and K^ DS ^ obtained with 
MC@NLO. This is not the case for the parton-level fixed-order computation, where 



8 In fact, Wt production is not implemented in HERWIG. We simply computed the LO matrix 
elements, and used the Les Houches interface [33] to give HERWIG the hard events, as in the case 
of MC@NLO. 
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the factors K display a significant dependence upon p^ e . This raises the following 
two issues. At small-p T , the fixed-order results lack Sudakov suppression at p\ — > 0; 
it follows that the predictions obtained for small values of p^ eto ' 1 are not particularly 
meaningful. At larger p T 's, the real-emission matrix elements dominate, and an NLO 
computation is expected to be reliable. Unfortunately, since the veto does not enter 
the LO results, it is not possible to sensibly estimate the impact of higher-order 
corrections: small changes to p^ eto ^ lead to large changes in the factors K, which 
may effectively mask the impact of tt interference. In summary, the imposition of a 
veto appears to be somewhat problematic in parton-level, fixed-order computations. 
This issue would imply a much (perhaps overly) larger theoretical systematic error, 
relative to an MC-based simulation, in the comparison between predictions and data. 

5.2 Gauge (in) dependence of DR cross sections 

As pointed out in sect. 4.2, the definition of the DR cross section, in which certain 
diagrams are removed from an amplitude, is not gauge invariant. In this section, we 
will argue that this is not a problem in practice, by repeating the DR calculation 
in a number of alternative gauges. We have checked several different observables, 
namely the total cross section; single-inclusive rapidity and transverse momentum 
distributions of the t, W and their decay products; and azimuthal and transverse 
momentum correlations between the charged leptons. We present here results for 
the total cross section and the transverse momentum p^ of the charged lepton pair. 
These are representative of the corresponding results for other observables. 

Our original calculation of the NLO cross section in the DR approach was carried 
out in the Feynman gauge, where the gluon propagator is given by: 

^ = (5-4) 

The sum over the polarisation degrees of freedom of initial-state gluons has been 
restricted to transverse polarisations. In order to check the sensitivity to the choice 
of gauge, we repeated the calculation in various covariant and non-covariant gauges. 
We summarise these results as follows. 

First, we note that the gauge dependence can only enter in the gg channel. This 
is because only diagrams involving two external gluons have broken gauge invariance, 
and these occur in the qg and gg channels. However, the qg channel has no resonant 
diagrams, thus no diagrams are removed. The diagrams kept in the gg channel 
(before crossing) are shown in figure 3. Only the third diagram in figure 3 leads to 
a gauge dependence, due to the presence of the gluon propagator. One can consider 
the general family of covariant gauges that result from replacing the propagator of 
eq. (5.4) with: 

^ = "^(^+(l-A)^). (5.5) 
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Figure 3: Diagrams kept in the gg channel after removing resonant contributions, shown 
before crossing. These formally represent also the two additional diagrams in which the 
two gluons arising from the fermion line have their momenta exchanged. 



This is useful, as A = 1 reproduces the Feynman gauge calculation. One thus has 
a continuous parameter in the DR calculation, such that it tends smoothly to the 
Feynman gauge calculation in a well-defined limit. However, it is easily seen from 
the form of the three gluon vertex that the second term in eq. (5.5) decouples from 
the amplitude. Thus, DR results do not depend on the choice of gauge within the 
family of covariant gauges. 

One may also consider non-covariant gauges, in which the gluon propagator is 
given by 9 : 



n ■ k (n ■ ky 



1. (5.6) 



The third term decouples for the same reason as in a covariant gauge, but a gauge 
dependence results from the second term. Here we report the numerical results for 
the MC@NLO cross section in DR, obtained with non-covariant gauges, for the three 
choices n 2 > 0, n 2 = 0, and n 2 < 0. We start in table 3 with the results obtained by 
integrating over the whole phase space, except perhaps for the veto imposed on the 
second-hardest B hadron of the event (see sect. 5.1). For reference, we also report the 
results obtained in the Feynman gauge (which is our default). We define a relative 
difference as follows: 

5 10 3 x <T<DR ' > ( non — covariant) — <j( DR ) (covariant) 

ct (dr) (covariant) 

We observe very small relative differences, which are actually compatible with zero 
within the statistical errors of our runs (each of which consists of 500k events). Note 
that this conclusion holds regardless of whether a veto is imposed or not. 

We have also checked that the same conclusion applies to more exclusive observ- 
ables. It is obviously impossible to reach the same level of statistical accuracy as 
that in table 3 for a differential distribution, especially in the tails of the transverse 



9 The gluon propagator also depends in general on a gauge- fixing parameter a, which we set 
equal to zero in what follows. See e.g. refs. [34, 35]. 
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Pt 


- 10 


J. vet °) 

FT 


= 30 


FT 


= 50 


FT 


= 70 


FT 


= oo 


a 


5 


a 


5 


a 


5 


a 


5 


a 


5 


Covariant 


34.66 





41.86 





44.61 





45.63 





46.33 





n 2 > 


34.70 


1.15 


41.90 


0.96 


44.70 


2.01 


45.73 


2.19 


46.42 


1.94 










n 2 -- 


= 










34.75 


2.60 


41.94 


1.91 


44.69 


1.79 


45.71 


1.75 


46.42 


1.94 


n 2 < 


34.70 


1.15 


41.87 


0.24 


44.61 





45.63 





46.37 


0.86 



Table 3: Gauge dependence of DR rates. The cross sections a are given in pb, and veto 
values p^ eto ^ in GeV. Note the factor 10 3 in the definition of S, eq. (5.7). The notation 
p(ve to ) _ qq denotes no veto at all. 



momentum spectra which fall rather steeply. We have computed the ratios 

= — T7^~ (non — covariant) / — ——(covariant) (5.8) 

bin-by-bin for the observables O (we have considered rapidities and transverse mo- 
menta, both for single-inclusive observables and for correlations), and found all ratios 
to be compatible with one. The typical value of \R(0) — 1 j is actually always much 
smaller than the statistical error affecting this quantity. The latter may be of the 
order of a few tens of percent in individual bins in the tails of p T distributions. 
However, since bin-by-bin fluctuations tend to integrate to zero if several contiguous 
bins are considered, we also fitted ratios in eq. (5.8) for those cases in which O is a 
transverse momentum, by assigning arbitrarily the same 0.1% relative errors to all 
bins, with the functional form 

oi + a 2 p T + a 3 p T . (5.9) 

This ensures that the tail of the p T distributions is treated in the fit on the same 
footing as the peak region, thus dealing in an efficient way with the problem of bin- 
by-bin fluctuations, and allowing us to uncover hints of non-flatness in R(p T ) as a 
function of p T . In this way, we have again found that all R(p T ) we have considered 
are compatible with one in the whole p T range. As an example, we present here the 
results obtained for the transverse momentum p^ of the lepton pair, for the four veto 
choices, and for the gauge n 2 > 0. The ratio R(p^) are presented in the upper pane 
of fig. 4. The result is also given for the case in which no veto at all is applied. As 
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DR(n 2 >0)/DR(covariant) 

—a a b n_ ..a^^a-^a a b a e a u ■ T3 — tl ff~ 
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Error on fit 



50 100 150 200 

p< u) (GeV) 

Figure 4: Gauge dependence of the DR cross section, as a function of the p T of the lepton 
pair, for a non-covariant gauge with n 2 > 0. Upper pane: results of the fits, using the 
form in eq. (5.9). Lower pane: envelope of curves obtained by varying all parameters of 
the fit within their error ranges (see the text for details), divided by the best fit curve, 
minus one. Black solid, red dashed, blue dotted, and green dot-dashed lines correspond to 
p(vet°) _ ^ anc | Q e y respectively. The magenta solid line with open boxes is 

obtained without imposing any veto. 

one can see from the figure, all curves resulting from the fit are remarkably flat (they 
cannot actually be easily distinguished). Some of them may be seen to depart from 
one at the largest values, but this behaviour is not statistically significant. To 
show this, we present in the lower pane of fig. 4 the envelope of the curves obtained 
by considering all combinations of values a>i = a° ± A(aj), with a° and A(a$) being 
the best values and their fitting errors respectively (as given by MINUIT [36]). The 
envelope is then divided by the best fit curve and unity is subtracted, so as to give 
an upper bound for the error affecting the fitting procedure. Only the result relevant 
to p^ eto ^ = 10 GeV is presented in the lower pane of fig. 4, since all the others are 
essentially identical. 

We therefore conclude that, regardless of the observable studied, the impact of 
gauge dependence in the DR computation can be safely neglected in the numerical 
studies that follow. 

5.3 Impact of interference 

In order to gauge how much of the difference between DS and DR is due to the 



1.2 



1.0 



0.8 



0.6 
0.10 
0.05 
0.00 
-0.05 
-0.10 



24 



interference term alone, it is useful to define the quantity: 

da™ = da® + A* - "J) #3 , (5.10) 

where the label NI stands for non-interference. Clearly 

da^ - d*™ = J2j , (5.11) 

which is then a direct estimate of the interference contribution, free of contaminations 
due to doubly- resonant contributions which are present in eq. (4.15). Although 
eq. (5.11) is not a physical quantity, it is useful to compute it because of the expression 
for da^-da^ in eq. (4.15). As can be seen there, this quantity is the interference 
term eq. (5.11), plus the difference between the doubly-resonant contribution and the 
subtraction term. Thus, eq. (5.11) would be identical to do~( DS ^ — d<j( DR ^ in the limiting 
case V = T> (which, we recall, is impossible to achieve if local cancellation and gauge 
invariance are simultaneously imposed). 

Results for the total cross sections after implementation in MC@NLO are shown 
in table 4, where we have defined the relative differences: 

ct (dr) ' V ) 

5 2 = r4 • (5.13) 

2 (j(DR) V 1 

Overall, we observe that the relative differences are not large, and decrease for tighter 
vetoes. This is to be expected given that by requiring a stricter veto, one filters out 
more of the tt process. Note that the relative differences due to the interference 
term alone are larger than those between DR and DS, implying that the terms X and 
T> — T> cancel out to some extent. 

Although not reported in this table, we have studied the scale dependence of 
these results, by varying the renormalisation and factorisation scales between m t /2 
and 2m t . The two scales are always set equal to a common value, which is somewhat 
restrictive but sufficient for our present purposes. Although the individual cross 
sections depend very mildly on scales (the variations w.r.t. the central values being 
-1%, ^2%, and tl% for DR, DS, and NI respectively), this is not the case for 
(j( DS ) _ (j( DR ) and <j( DS ) — &( m \ which have variations w.r.t. their central values of 
about ±20% and ±30% respectively. This behaviour is not surprising: the individual 
cross sections have been defined with the specific purpose of describing the Wt cross 
section, which is known from the literature to have a mild scale dependence. On 
the other hand, as explicitly shown in eqs. (4.15) and (5.11), the differences of cross 
sections considered in table 4 are dominated by LO tt production, which has a 
significant scale dependence (see e.g. ref. [37] for a recent update). 
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ct (dr) 


(j(DS) 




(j(DS) _ a (DR) 


Si 


a (DS) _ a (m) 


6 2 


p{ veto) = 1Q QeV 


34.66 


33.89 


35.08 


-0.77 


-2.2% 


-1.19 


-3.4% 


= 30 GeV 


41.86 


40.74 


42.83 


-1.12 


-2.7% 


-2.09 


-5.0% 








(veto) - 

p y T = 5 


GeV 






44.61 


42.92 


46.59 


-1.69 


-3.8% 


-3.67 


-8.2% 


p[ veto) = 7Q 


45.63 


43.65 


48.24 


-1.98 


-4.3% 


-4.59 


-10.0% 








FT - 


= oo 






46.33 


44.12 


49.58 


-2.21 


-4.8% 


-5.46 


-11.8% 



Table 4: Results for the total DR, DS, and NI cross sections (in pb) denned in eqs. (4.9), 
(4.14), and (5.10), for various values of the veto. We also give the relative differences as 
defined in eqs. (5.12) and (5.13). The notation p^ eto - ) = 00 denotes no veto at all. 



We have also carried out similar studies on differential distributions. The results 
presented in table 4 are, not surprisingly, dominated by the small transverse momen- 
tum regions, where the bulk of the cross section lies. Analogously to what was done 
in eq. (5.8), we have computed the ratios 

, \ dcr (DS) /r/rx (DR) , s r/rr (Nl) / 

-icr / -icr • ^"(^^r/V- <" 4 > 

for several observables. In the cases in which O is a rapidity, these ratios are flat 
over the whole kinematically-accessible range, with values consistent with those given 
in table 4. The situation is more interesting if one considers transverse momenta, 
since in such cases the ratios in eq. (5.14) display a non-trivial shape. In particular, 
the largest values of | — 1 1 and | — 1 1 , amongst the observables we have 
studied, are found in the tail of p^ (which is the reason why this observable has 
been considered as a case study in this paper). 

We start by presenting results for the quantities defined in eq. (5.14). In order 
to be able to superimpose five curves on the same plot and keep visibility (since 
histograms would blur the picture because of bin- by-bin fluctuations in the tail), 
the ratios have again been fitted with the functional form of eq. (5.9). The results 
are given in the upper panes of fig. 5. As can be seen from the figure, the impact 
of interference can be very large for large enough values of p? l \ regardless of the 
choice for p^ eto ^ (although, of course, it is less significant for small vetos). It must be 
stressed however that, at large p^ l \ the cross section is small. This is documented 
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Figure 5: Upper panes: results for the ratios denned in eq. (5.14), as a function of and 
for various vetos. Lower panes: relative scale dependence (see eq. (5.16)). The linestyles 
are the same as those of fig. 4. 




pj""°> = 10 GeV 




Figure 6: Left pane: differential DR, DS, and LO distributions in Px . Right pane: 



integral of the same distributions in the range pi£ nm * > < p^ l> < oo, divided by the respective 



total rates. These results are relevant to the case p T 



(veto) 



10 GeV. 



in fig. 6, where we present in the left pane the differential distributions in Px , as 
computed with DR, DS, and at LO. The same information is presented in the right 
pane of the figure, in an integral form: 

J-Eb*-') = J- /" tf'-^j. (5.15) 

The results of fig. 6 have been obtained by choosing pi! eto ^ = 10 GeV. Although 
the absolute value of the differential cross section has a non-negligible dependence 
on p^ eto \ its shape is relatively stable against variations of p^ eto \ Thus, the right 
pane of fig. 6 can be used with table 4 for estimating the number of events with 
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Pt > Pt for any choice of the veto (or in its absence). 

The comparison between fig. 5 and fig. 6 shows that, for the majority of events, 
the impact of interference is moderate. On the other hand, these results imply that 
conclusions concerning interference of Wt with tt are observable dependent. One 
particular analysis may be sensitive to an observable which receives much larger 
contributions from J and V — V than the total rate. With fig. 5 we have presented 
a worst-case scenario among the observables studied. 

It is also interesting to study the scale dependence of these results. This cannot 
be done in the same way as for the total rates of table 4. There, differences of cross 
sections have been considered, and the relative variation w.r.t. the central values 
computed. In the case of the tail of a steeply falling differential distribution, bin-by- 
bin differences lead to such small numbers that it is impossible, in practice, to obtain 
a statistically significant result for the relative variation, since the latter would be a 
ratio of two extremely small numbers. On the other hand, we are able to compute 
the scale dependence of the ratios defined in eq. (5.14). In the lower panes of fig. 5, 
we present the quantities 



(and the analogous one for NI) as a function of . Equation (5.16) will not give 
directly the scale dependence of a tt-dominated cross section, as was the case for 
the differences of rates of table 4, but rather that of the fraction of the DS and NI 
cross sections arising from tt-like contributions. As can be seen from the figure, the 
scale dependence of such a fraction is largest in the region where the interference 
between singly- and doubly- resonant diagrams is largest. This is consistent with the 
naive expectation that I^t-like contributions have a much milder scale dependence 
than tt-like ones. On the other hand, the study of scale dependence shows that the 
qualitative conclusion that one draws from fig. 5 is perturbatively stable, and that 
the interference is significant only in the large-p T region. 

6. Discussion 

In this paper we have presented the first calculation of single-top production in 
association with a W boson, accurate to NLO in QCD and interfaced with parton 
showers according to the MC@NLO formalism. In studying the Wt channel, our 
main aim has been to investigate whether or not one can define it as a separate 
production process, since only in the case of a positive answer is the computation 
of the full NLO corrections meaningful and feasible with present technologies. The 
NLO+parton shower framework provides an excellent platform for addressing the 




-1, 



k = 1/2,2 



(5.16) 
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problem. Parton showers give the capability of studying final states in a realistic 
environment, and testing ideas that can be applied without any modifications to 
an experimental analysis. The underlying NLO computation is, on the other hand, 
directly sensitive to interference effects with tt production. 

We have considered two separation mechanisms for defining the Wt channel: 
diagram removal (DR) and diagram subtraction (DS). We have defined them in such 
a way that it is not necessary, in order to generate full hadronic events with high 
efficiency, to apply any kinematic cuts. In this way complete flexibility is achieved, 
to study the effects of any cuts whose aim is that of enhancing the Wt "signal" w.r.t 
the tt "background". As an example, we have considered only one such cut, namely 
a veto on the transverse momentum of the second-hardest B hadron of the event. 

The DR cross section is obtained by removing the doubly-resonant diagrams from 
the calculation of the Wt amplitude, and thus violates gauge invariance. We have 
demonstrated that this is not a problem in practice. The DS cross section involves 
a modification of the naive Wt cross-section by a gauge invariant subtraction term 
which acts to remove the doubly-resonant contribution from tt final states. Given 
that this is done at the cross-section level, an interference term between the tt and 
Wt processes remains. 

Ultimately, we are in a position to decide whether it is possible, and useful for 
data analysis, to define the Wt channel as a production process in its own right. The 
key question is whether an approximate treatment of interference within a higher- 
order computation has to be preferred, for the sake of a counting experiment, to 
an exact treatment of interference without QCD corrections. Furthermore, the sep- 
aration of Wt and tt processes allows one to consider NLO Wt results (which are 
0(9w a D) alongside NLO tt results (which are Ofaf)). This makes a realistic descrip- 
tion of phenomenology possible (since it is well known that higher-order corrections 
for tt production are crucial in this respect), in spite of being unable to treat the 
problem exactly in its full complexity. In this paper we have found that NLO cor- 
rections are larger than interference effects, especially if cuts are employed that are 
designed to enhance the Wt signature. The answer to the above question there- 
fore appears to be yes. However, the definitions that render it possible to compute 
higher-order corrections also imply an ambiguity in the theoretical predictions, that 
we identify with the difference between DR and DS cross sections. We have found 
that such an ambiguity is observable-dependent, and in particular is larger in those 
regions of the phase space that correspond to large transverse momenta. Since this 
ambiguity will be directly related to the theoretical systematic error, it follows that 
the accuracy with which Wt properties can be measured depends on the observables 
that are most relevant to a given analysis. 

To conclude, having implemented the Wt production mode in an NLO plus 
parton shower context, we find based on our subsequent analysis that it does indeed 
seem feasible to analyse this process at the LHC. 
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A. W~t and W + i production 

In this appendix we demonstrate that the squared amplitudes for single top produc- 
tion with an associated W boson are independent of whether the final state top is a 
quark or an antiquark. 

Firstly there is the issue of a possible charge asymmetry due to the electroweak 
coupling of the W boson. This merely constrains the helicity of the fermion line 
which includes the top quark in the graph, which occurs either as an open line or a 
closed loop. The t — > t transformation changes the direction of the line, and leads to 
a constraint on the helicity of the corresponding antifermion line. However, the sign 
and size of the coupling is the same in both cases. 

Furthermore, there is the issue of a possible QCD asymmetry. One must consider 
in more detail the open or closed fermion loop associated with the top quark (and 
the b quark from top decay). The replacement t — > i can affect the amplitude in two 
ways: 

1. It changes the sign of all momenta which occur in uncut propagators along the 
top quark line 10 . If Pi are the propagating momenta, the fermion trace in each 
diagram will contain the following terms: 

{i> x -m t )(tf 2 -m t ) ■■■{i>M -mt), (A- 1 ) 

where M is the number of uncut fermion momenta, and there may be additional 
Dirac matrices between the propagator factors. Interchanging top and antitop 
quarks gives instead the terms 

(- i> x -m t )(- i> 2 -m t ) jf M -mt) = (~ l ) M (^i+m t )(^2+m t ) ■ ■ ■ (0M+m t )- 
(A.2) 

10 Cut propagators correspond to terms in the fermion trace of form $±m t depending on whether 
quark or antiquark spinors are involved i.e. the sign of the momentum is not changed. 
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Given that the number of Dirac matrices in each interference graph is even, 
only terms involving an even number of masses survive. Then the amplitude 
for antiquark production is related to that for quark production via 




(A.3) 



where At,t is the amplitude associated with a given interference diagram. Here 
Cdiag is the parity under t — > t of the diagram i.e. Cdiag = ±1. 

2. The colour factor may be affected. In each of the diagrams for Wt production, 
one always has 



where C 4i j are the colour factors for the amplitude with a top and antitop quark 
respectively, and C co \ — ±1. 

The total parity under t — > t is then given by 



Squared real emission diagrams will always have an even number of uncut prop- 
agators, and also symmetric colour factors under t — > t. Regarding interference 
diagrams, there are 48 in total. These can be subdivided into gg, bb, bg and bb initial 
states. Then diagrams with other quarks in the initial state form a subset of those 
already specified. The gg and bb diagrams are always associated with real emissions. 
For the bb (and, hence, qq) initial state, there is only one Feynman amplitude and 
hence no interference term is possible. The qg initial states can be associated with 
real or virtual emissions. By evaluating the number of uncut fermion propagators 
and the colour factor for each graph, one can find its parity under t — > t using 
eq. (A. 5). 

There are two types of diagrams: 

1. Diagrams with no triple gluon coupling. These all have symmetric colour fac- 
tors under t — > t, and an even number of uncut fermion propagators. Hence 
C — 1 for such graphs. 

2. Diagrams with a triple gluon coupling. These all have antisymmetric colour 
factors, and an odd number of uncut fermion propagators (this latter fact can 
be easily appreciated by considering removing a gluon line from a fermion in 
graphs of type 1 and reattaching it to a gluon line). Hence C = (— l) 2 = +1 
for these graphs. 



Ct — C co iCt, 



(A.4) 



C — CdiagCcol — ±1- 



(A.5) 
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One finds that every graph is even upon replacing top quarks by antitop quarks, and 
so the total squared amplitude for Wt production is the same for both t and t. 

Even if W~t and W + i matrix elements are identical at this order, there is still 
the possibility of a charge asymmetry arising from differences between the b and b 
parton densities 11 . However, these densities are equal in the global fits available at 
present. Thus, overall, the cross- sect ions and differential distributions considered 
here for top and anti-top production are strictly equal. 



B. Calculation of A4 in the helicity formalism 

NLO subtraction formalisms such as FKS achieve a good numerical stability through 
the local cancellation between large contributions of opposite sign, which correspond 
to the real-emission matrix elements in the soft and collinear regions, and the coun- 
terterms. The counterterms are in fact constructed by calculating the above limits 
of the matrix elements, whose form is universal. If one considers e.g. the initial-state 
collinear limit p 2 \\ k of the process 

a(pi)+/3(p 2 ) — X(K)+S(k), (B.l) 

with X a set of final state particles which are not relevant here (in our case, X = 
t + W), the matrix element squared is (see e.g. eq. (B.41) of ref. [15]) 

M( Pl ,p 2 ) ^ ^ \p(z)M^( Pl ,z P2 )+Q(z)M( Pl ,z P2 )] , (B.2) 
k- P2 i J 

where P are the usual Altarelli-Parisi kernels, Q are other universal kernels (given 
at the leading order in eqs. (B.42)-(B.45) of ref. [15] for initial-state collinear split- 
tings, and in eqs. (B.31)-(B.34) of that paper for final-state collinear splittings; these 
kernels are thus different for spacelike and timelike branchings already at the lead- 
ing order), Ai^ is the relevant Born contribution, and Ai is a Born-like function, 
which however keeps track of the azimuthal correlations in the branching process. 
The contribution of QAi vanishes upon integration over the azimuthal angle of the 
branching, which is why this term can be neglected in the analytical computation of 
the collinear divergences in 4 — 2e dimensions. Locally, it is different from zero if the 
parton involved in the branching which enters the hard reaction is a gluon, and there- 
fore needs to be taken into account for the construction of the local counterterms in 
a (efficient) numerical NLO program. 

Unlike the case of s- and t-channel single-top production, the Wt channel has a 
gluon entering the hard reaction at the Born level, and we thus need to compute the 
relevant At (which is process dependent, whereas the kernels Q are universal). This 



11 As already discussed in refs. [7, 11]. 
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+ 





Figure 7: Diagrams used in the calculation of M for — > +5(^2) + (P3) +*(P4)j 
i.e. all momenta are denned to be outgoing. 



term is (see eq. (B.23) of ref. [15]) 



M = FRe 




(B.3) 



with A± the Born amplitude for Wt production with a positive/negative helicity 
initial state gluon respectively, and T a term which includes the flux and spin- and 
colour- average factors, and the coupling constants. We have computed the r.h.s. of 
eq. (B.3) in two independent ways, and found agreement. In this appendix we present 
some details on how the computation can be carried out by using spinor helicity 
methods [38, 39, 40, 41, 42]. As is conventional when applying such methods, we 
take all momenta to be outgoing in this appendix, in contrast to the rest of the paper. 
Furthermore, we relabel the outgoing momenta k 2 , k\ as p 3 and p 4 respectively (see 
fig. 7), in order to be able to use a compact notation in what follows. 

The relevant diagrams are illustrated in fig. 7. Besides the null momenta p\ and 
p 2 , we define the null vector 



where U\ = u — m 2 = {p 2 + pi) 2 — m\. As is common, we adopt the notation 



2 

mi 

P5=PA P2, 

Itl 




A; ± )=m ± (A;), (k 



±1 — 



u±(k) 




for massless quark spinors. Defining the Mandelstam invariants 



s = (p 1 +p 2 ) 2 - t = (p 1 +p 4 ) 2 , 



(B.6) 
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the contributions A^ and of the positive helicity diagrams (denoted by (a) and 
(b) respectively in fig. 7) are given by (neglecting colour factors) 



(s) 2iV2 



A + = -j^fMMp*hT-)(^)(2+Wp&, (b.7) 

A + = ^j eMU{P4) [|l+>[25]<5-| +m|2->(l-|]y«(pi), (B.8) 

where e^ps) is the polarisation vector of the W boson, and we have used the standard 
notation 

[ij] = (i + \j-), (ij) = (i-\J+)- (B.9) 
In deriving eqs. (B.7) and (B.8) we have used the following expressions for the gluon 
polarisation vectors 

e + AP^)- ^ {q2) , e.j.M- . (B.10) 

Here q is an arbitrary null reference momentum, which has been set to p 1 . One is 
in principle able to choose a different reference momentum for the negative helicity 
diagrams, owing to the fact that helicity amplitudes are separately gauge invariant. 
However, if one also chooses q2 = P\ in the negative helicity case then diagram (c) 
in fig. 7 vanishes. This relies on the helicity of the light quark line being fixed by 
the chiral boson coupling. Then the negative helicity amplitude is given solely by 
diagram (d): 



11 

— 12+)[12](2-| + |2+)[15](5-| +m t |l-)(2-| 



y*u(pi). 



(B.H) 

It can be checked that \A^\ 2 + |^L 0) | 2 , after summing over the W boson and quark 
spins, is equal to the known Born result, eq. (3.5). For the helicity interference term, 
one uses eqs. (B.7), (B.8), and (B.ll) to find (A^ = A { + } + A { + ] ): 

= 8 [l %^\ ml - m 2 t )(2m> w + ml). (B.12) 

Note the presence of squared spinor products, which cannot immediately be evaluated 
to form dot products. This is to be expected, given that M, is not a Lorentz invariant 
quantity. To evaluate the quantities in eq. (B.12) one can parameterise 4- momenta 
in the center of mass frame of the incoming particles as follows: 



/ (^,0,0,^): (B.13) 



P2=( #0,0,-^1; (B.14) 



2 1 1 1 2 

p 5 = (E,Esin8cos<f),Esm6sm4>,Ecos8) . (B.15) 
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Basis spinors satisfying the massless Dirac equation with conventional normalisation 
vf(pi)u(pi) = 2p® are (choosing arbitrary phases) 








\ 



-y/E(l - cos6)e~ i4 
y ,/E(l + cos6) ) 



u\{pi) 



( s*5 x+ \ 





U\{P2) 



( \ 

s*5x+ 



(B.16) 



/ 



from which one finds 

(kp 2 ) [15] 2 (25) 2 



Re 



[kp 2 ] [12]2 



(15)(25) 



[2 cos 2 (0-0 2fe ) - 1], 



(B.17) 



with 02fc the azimuthal relative angle between k and p 2 . In the degenerate kinematics 
of the collinear limit, we may choose to absorb in a redefinition of <f) 2 k- Combining 
these results with eq. (B.12), and reinstating the colour, flux, average factors and 
coupling constants, gives finally 



M 



16iVs 2 u\rriw 



(2cos 2 02fc - l)(ravf _ m t)( 2m w + m t)( m w m t - ut )- (B.18) 



Note that this does indeed give zero contribution to the cross-section when integrated 
over the full domain of the azimuthal angle 2 fc relevant to the branching, as discussed 
above. 
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